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fSj . Abstract 

^-y We have studied the famous classical pseudotensors in the small region limit, 

both inside matter and in vacuum. A recent work [Deser et al. 1999 CQG 16, 
2815] had found one combination of the Einstein and Landau-Lifsfntz expres- 
sions which yields the Bel-Robinson tensor in vacuum. Using similar methods 
we found another independent combination of the Bergmann- Thomson, Papa- 
petrou and Weinberg pseudotensors with the same desired property. Moreover 
£> ■ we have constructed an infinite number of additional new holonomic pseu- 

dotensors satisfying this important positive energy requirement, all seem quite 
artificial. On the other hand we found that Moller's 1961 tetrad-teleparallel 
I/"") ■ energy-momentum expression naturally has this Bel-Robinson property. 
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1 Introduction: quasilocal quantities for small re- 
gions 



The localization of gravitational energy-momentum remains an important problem in 
GR. Using standard methods many famous researchers each found their own expres- 
sion. None of these expressions is covariant, they are all reference frame dependent 
(referred to as "pseudotensors"). This feature can be understood in terms of the 
equivalence principle: gravity cannot be detected at a point, so it cannot have a 
point-wise defined energy-momentum density. Now there is another way to address 
the difficulty 

The new idea is quasilocal: energy-momentum is associated with a closed 2 surface 
surrounding a region pQ. A good quasilocal approach is in terms of the Hamiltonian 
Then the Hamiltonian boundary term determines the quasilocal quantities. In 
fact this approach includes all the traditional pseudotensors jSlHj- They are each gen- 
erated by a superpotential which can serve as special type of Hamiltonian boundary 
term. 

A good energy-momentum expression for gravitating systems should satisfy a va- 
riety of requirements, including giving the standard values for the total quantities for 
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asymptotically flat space and reducing to the material energy-momentum in the ap- 
propriate limit. No entirely satisfactory expression has yet been identified. One of the 
most restrictive requirements is positivity. A general positivity proof is very difficult. 
One limit that is not so difficult is the small region limit. This has not previously been 
systematically investigated for all the classical pseudotensors expressions. In matter 
the expression should be dominated by the material energy-momentum tensor. In 
vacuum we require that its Taylor series expansion in Riemann normal coordinates 
have at the second order a positive multiple of the Bel-Robinson tensor. This will 
guarantee positive energy in small vacuum regions. 

2 The classical pseudotensors 

Recall the Einstein field equation G^ u = R^ u — ^g^R = kT^ u . Matter energy- 
momentum has a vanishing covariant divergence V^T^j, = 0, but in curved spacetime 
this is not in the form of a conserved energy-momentum relation. But one can rewrite 
it in the form of a divergence 

d^g{T\ + t\) = ^ (1) 

here t^ v is the gravitational energy-momentum pseudotensor. The energy-momentum 
complex T M j, is then given as 

kT» u = K^g~(T» u + t\) = d x U^ x \ (2) 

where Uj 1 ^ is called the superpotential. 

3 Riemann normal coordinates and the adapted 
tetrads 

For the total quantities of an isolated gravitating system the various expressions give 
the expected weak field asymptotic values. However they are quite different in the 
strong field region. To study the quasilocal quantities one can Taylor expand the 
Hamiltonian, including the divergence of its boundary term, in a small spatial region 
surrounding a point. The reference is the flat space geometry at this origin. Riemann 
normal coordinates (RNC) satisfy 

fiV(°) = V^u, <9 A sv(0) = 0, d 2 flu g a(3 (0) = --(R a n/3v + #cm^)(0). (3) 
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Later we will also need the associated adapted orthonormal frame (aka tetrad, vier- 
bein) which satisfies g^ u = r] ab e a M e fe v and 

e%(0) = ^, d v e a v {0) = 0, I%(0) = 0, 3^(0) = ^VC )- ( 4 ) 

For the energy-momentum density we want at non-vacuum points the results to be 
dominated by the material energy-momentum. At vacuum points one wants a result 
proportional to the Bel-Robinson tensor. That will guarantee the proper positive 
energy property; i.e. the associated energy-momentum vector will then be future 
pointing and non-space like. 



4 Quadratic curvature combinations 

The Bel- Robinson tensor and the tensors S a p^ V) K a p^ v and T a ^^ v are defined 

as follows 
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A recent work [3] had found exactly one pseudotensor expression, a certain combina- 
tion of the Einstein and Landau-Lifshitz expressions, which yields the Bel-Robinson 
tensor in vacuum. They argued that this combination is unique under their assump- 
tions. 



5 Classical holonomic pseudotensors 

The well known classical superpotentials associated with the Einstein, Landau-Lifshitz, 
Bergmann- Thomson, Papapetrou, Weinberg and M0ller(1958) energy-momentum com- 
plexes are 

(9) 

-g L U^\ (10) 

(11) 



pi? HM = V=gg ma g nd €d 
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= (v m c V nd - \v mn V cd ) 9c d V ah &nl (12) 

58 ^ [H = ^(r^-r^). (13) 

The pseudotensors are obtained according to one of the prescriptions: 

T a » = d u Uj H , T aix = d v U a[lxv \ = d 2 al3 H^ v ®. (14) 

Inside matter at the origin, the RNC expansion results are 

£/(0) = 5/(0) = P/(0) = W/(0) = 2G Q » = 2kT/(0). (15) 

This is as expected from the equivalence principle. However the M0ller(1958) expres- 
sion gives 

M/(0) = Rj(0) = k (Tj - 1 V Jt) (0) + 2«T/(0). (16) 
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This result is not acceptable. From [Sj, the small vacuum {G^ u = 0) region non- 
vanishing terms are 

Ej = -2T Xaa T^ + 6^T XaT T^, (17) 
L at3 = Y a Xa (T l3X ' 7 -T Xal3 )-T X a a (T pX(T + T rT ^)+g a ^T XaT T aXT . (18) 

Using a similar technique we found, for the other pseudotensors (note Papapetrou 
and Weinberg depend explicitly on the Minkowski metric r] a p), 

pap = L aP + h Xa( r aP x ^ + r Pa x j_ h Xf3^ ra c, x + r ** x ^ ^ 

W af3 = -2r Aa Q r A ^ + ^r ACTr r A ^-^V^^ (20) 

Mj = r ACTa r A ^-r ACTa r^ + (/V A -/V a )^,A.. (21) 

Here h a p := g a/ 3 — r] a p and h a p = —^R a x/3aX X x u + 0(x 3 ) in RNC. From [5] in vacuum 
we have 



9l u E Q/3 (0) - -(AB a)3flu - S a/3fMU ), (9^L Q/3 (0) - - (7B a/3flu + -S a p^ u ) . (22) 



In order to obtain the Bel- Robinson tensor, Deser et al. used a "by hand" combination 

^Iv (nE a /3 + L a f^\ = B^^, (23) 
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Using similar methods we obtained for the other pseudotensors in vacuum at the 
origin 
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(25) 
(26) 



(Here we have included for completeness the result of the M0ller(1958) expression, 
even though it does not have a good matter interior limit.) From this we find one more 
independent combination of the Landau-Lifshitz, Papapetrou and Weinberg pseu- 
dotensors with the same desired Bel-Robinson property. (Note: the earlier work cited 
above did not get this result, as they had excluded the explicit use of the Minkowski 
metric in the superpotentials they considered.) Inside matter and in vacuum at the 
origin, respectively, we find 



2£ Q /3 + -(3P Q/ g 
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2G a p(0), —d n U 1L a p + -(3P Q/ 3 
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W, 



a/3 J 



(0) — B a [3^ u . 

(27) 



6 A large class of new pseudotensors 

Moreover we have constructed an infinite number (a 3 parameter set) of new holo- 
nomic pseudotensors all satisfying this important Bel-Robinson/positivity property. 
The new general superpotential is 

7/ [H _ TT [H _i_ J c ih^^a U -n + C2h tJl7T T u an + Czh^T^a 1 ,no\ 

where C\ to C5 are constants and h a p := g a p — r} a p. In Riemann normal coordinates 
h a P = \R a ^ K x^x K + 0(x 3 ). Actually, the leading superpotential E U a [llu] in flSE) is 
not necessary Freud's, it can be replaced by any other which offers a good spatially 
asymptotic and small region material limit. The resultant energy density inside mat- 
ter at the origin is, as expected, 

2re£/(0) = 2G/(0) = 2kT/(0). (29) 
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The RNC second derivatives in vacuum at the origin are 

(4 - 2ci + c 2 + c 3 - 4c 4 + 3c 5 )B a/3 ^ v 
dL£ af) (0) = J I -\{2 - Cl - 4c 2 + 5c 3 - 2c 4 - 3c 5 )S, 



a/3/j.u 



+(ci + c 2 - 2c 3 - 2c 4 )K, 



a/3[ii/ 



(30) 



where c 4 = c 4 — |c 4 . Note that we can choose 



4-2ci + C2 + c 3 -4c 4 + 3c 5 >0, (31) 
2 - ci - 4c 2 + 5c 3 - 2c 4 - 3c 5 = 0, (32) 
ci + c 2 - 2c 3 - 2c 4 = 0. (33) 

The solutions can be parameterized as follows, 

ci+c 2 -2c 3 < 1, (ci, c 2 , c 3 , c 4 , c 5 ) = ^ci, c 2 , c 3 , i(ci + c 2 + 2c 3 ), ^(2 - 2c 4 - 5c 2 + 7c 3 ^J . 

(34) 

Then the second derivatives of the new pseudotensors in vacuum are 

3LA/?(0) = ^(1 - Cl - c 2 + 2c 3 )B af3 ^. (35) 

As there are three arbitrary constants that one can tune, obviously there exists an 
infinite number of solutions with any positive magnitude of B a p^ v . They all appear 
to be highly artificial. It seems that there is no obstruction in going to higher order. 
From our analysis we infer that there are an infinite number of holonomic gravitational 
energy-momentum pseudotensor expressions which satisfy the highly desired small 
region Bel- Robinson/positive energy property. 



7 M0ller's 1961 tetrad-teleparallel energy- momentum 
tensor 

On the other hand, unlike the aforementioned mathematically and physically con- 
trived expressions, we found that M0ller's 1961 teleparallel-tetrad energy-momentum 
expression naturally has the desired Bel-Robinson property. The superpotential has 
the same form as Freud's but the indices now refer to a tetrad: 

e iU a [bc] = V=gg df r l fe s<;z- (36) 
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Expressed in terms of differential forms we have 

m a b Vb = d(T b c A Vb c a ). (37) 

The RNC and adapted frame expansion results inside matter, and the vacuum second 
derivatives at the origin, respectively, are 

m a b (0) = 2G a \0) = 2/tT a fe (0), d%m ab {Q) = -B^. (38) 

Thus the desired Bel-Robinson property is naturally satisfied. An important conse- 
quence is that the gravitational energy according to this measure is positive, at least 
to this order. (We expected this positivity result since in fact M0ller's 1961 expression 
has an associated positive energy proof (BJ.) Once again M0ller's 1961 tensor stands 
out as one of the best descriptions for gravitational energy-momentum. 
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